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Abstract 

We apply the eigenvectors from a variational analysis in lattice QCD to successfully extract the wave function of 
the Roper state, and the next Pn state of the nucleon, associated with the A^*(1710). We use the 2 + 1 flavour 
32^ X 64 PACS-CS configurations at a near physical pion mass of 156 MeV. We find that both states exhibit a structure 
consistent with a constituent quark model. The Roper d-quark wave function contains a single node consistent with 
a 25 state, and the A^*(1710) contains two, consistent with a 35 state. A detailed comparison with constituent quark 
model wave functions is carried out, obtained from a Coulomb + ramp potential. These results validate the approach 
of accessing these states by constructing a variational basis composed of different levels of fermion source and sink 
smearing. Furthermore, significant finite volume eff'ects are apparent for these excited states which mix with multi- 
particle states, driving their masses away from physical values and enabling the extraction of resonance parameters 
from lattice QCD simulations. 
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1. Introduction 

The wave function (and probability distribution) of a 
particle in a potential are fundamental to the very na- 
ture of quantum mechanics. In the non-re lativistic case, 
the entire spectrum of the particle can be determined by 
solving the Schrodinger equation. In quantum field the- 
ory, a Schrodinger-like probability distribution can be 
constructed for bound states by analogy, in the form of 
the Bethe-Salpeter amplitude [1 1. 

By considering a non-relativistic constituent quark 
model, the probability distributions of quarks within 
hadrons can be determined using a one gluon exchange 
potential augmented with a confining form ||2l [31. 
These models have been the cornerstone of intuition 
of hadronic probability distributions for many decades, 
and have been complemented with features such as 
meson-cloud dressing. Quark probability distributions 
can be calculated directly (without modeling) via Lat- 
tice Gauge Theory ID. VisuaUzations of the probabil- 
ity distribution on the lattice have been used to observe 
interesting physical effects such as Lorentz contraction 
lUllSl, quarks aligning with a magnetic field and diquark 
clustering |7|. Furthermore, the probability distribution 
can be used as a diagnostic tool, allowing finite volume 
effects and other lattice artifacts to be easily visualized 



and understood fSl . 

Hadron spectroscopy is a highly complex problem. 
Though it is relatively simple to see higher energy reso- 
nances of hadrons in colliders, apart from simple quan- 
tum numbers, properties more fundamental to the nature 
of these resonances remain elusive to experiment. 

Robust methods have been developed that allow the 
isolation and study of states associated with these res- 
onances in Lattice QCD HIlOliniiniinilllllSlIEl 
[iTlIIHl. In this study, we apply the variational method 
lfT9l l20l to extract the ground state and first two Pi i 
excited states of the proton associated with the Roper 
Ell and the A^*( 1710). We then combine this with lat- 
tice wave function techniques to calculate the probabil- 
ity distributions of these states at near-physical quark 
masses. We use the 2 + 1 flavour 32"* x 64 PACS-CS 
configurations ll22l at a near physical pion mass of 156 
MeV. 

The wave function of a hadron is proportional to the 
parity-projected |23| two-point Green's function, 

Gfjip, t)^Yj '''''^ (^" ± 1> <"l TiXii^XjiO) } |fi), 

(1) 
where Xi are the hadronic interpolating fields. In the 
case of the proton the most commonly used interpolator 
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is given by 



Xi(x) 



' (ul(x) Cys dt(x)) Uc(x), 



with the corresponding adjoint given by 

Xi(0) = e""^ ( db(0) Cj, mJ(0) ) «,(0), 



(2) 



(3) 



where u and d represent the up and down quark fields 
respectively. In order to construct the wave function, the 
quark fields in the annihilation operator are each given 
a spatial dependence, 

X\ (4 y, t -w) = e"*" {u\{x+y) Cys dh{x + z)) ujx + w), 

(4) 
while the creation operator remains local. This gen- 
eralizes G(p, f) to a wave function proportional to 
G{p,t;y,z,w). In principle, we could allow each of 
these coordinates, y, z, w, to vary across the entire lat- 
tice, however, we can reduce the complexity by tak- 
ing advantage of the hyper-cubic rotational and trans- 
lational symmetries of the lattice and considering the 
system centre of mass. A description of the probability 
distribution of a particular quark within the proton can 
be formed by holding the spatial location of two of the 
quarks fixed and calculating the third quark's amplitude 
at every lattice site. We focus on the probability distri- 
bution of the d quark from Eq. (pi, with the w-quarks 
fixed at the origin, i.e. 

Xi(x, 0, 1 0) = e"'' ( uli^ Cy, dh(^ + f) ) M,(f). (5) 

The coordinate x then represents the centre of the sys- 
tem, providing an origin analogous to a potential well 
in the non-relativistic case. The resulting construct is 
gauge dependent so we choose to fix to the Landau 
gauge. 

The variational method 1191 l20l is a well-established 
method [24] for extracting the excited state spectra of 
hadrons. Noting that the only time dependence in the 
two-point function lies in the exponential, we are able 
to construct the following relation, 

Gijdo + At) u"j = e-'""*' Gijito) u'] , (6) 

where m" are the right eigenvectors of the eigenvalue 
equation, 

( G-\to) G(to + At) )ij mJ = e^"'"^' Mf . (7) 

Similarly, we can construct the left eigenvector equa- 
tion. 
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Figure 1 : Tlie effective mass of tiie ground and first two excited states 
of the nucleon projected with only the right eigenvector. Effective 
projection is accomplished up to 3 time slices after the source, and the 
masses match those calculated in Ref . 1 22 1 . 



To project a single state we apply the eigenvectors to the 
parity-projected variational matrix. 



vfGtj(t)i/j^6a,3e 



-mt At 



(9) 



yf(G(fo + AOG-'(fo)),7 = e 






(8) 



The eff'ective mass can then be calculated 
from the projected two-point functions as 
m(t) = log(G(0/G(f +1)). While the eff^ective 
mass is insensitive to a wide range of variational 
parameters ifTTl . we follow Ref. ifTTl and select fo to be 
2 time slices after the source with At - 2. 

Different interpolators exhibit different couplings to 
the proton ground and excited states and hence can be 
used to construct a variational basis. The limited num- 
ber of local interpolators restricts the size of the op- 
erator basis ID. To remedy this, one can exploit the 
smearing dependence of the coupling of states to one or 
more standard interpolating operators in order to con- 
struct a larger variational basis where the^, and;^') from 
Eq. (fTb contain a smearing dependence. This method 
has been shown to allow access to states associated with 
resonances such as the Roper, A^*(1710) |17| and the 
A(1405) Il25l. 

The non-local sink operator used to construct the 
wave function is unable to be smeared, such that the 
standard technique of Eq. (J9]l cannot be applied. How- 
ever, Eq. (J6]l illustrates it is sufficient to isolate the state 
at the source using the right eigenvector Thus, the prob- 
ability distributions are calculated with each smeared 
source operator and the right eigenvectors calculated 
from the standard variational analysis are then applied 
in order to extract the individual states. As demon- 
strated in Fig.[T] clean projection of the first two excited 
states is obtained. We note how the plateaus commence 



at t = to - 2, where the correlation matrix analysis has 
been applied. 

In summary, the wave function for state a observed 
at Euclidean time t is 

X 

tr (ro ± 1) <Q| T{xi(xJ,t 0)xj(O) } |Q> mJ . 

By averaging over the equally weighted {U} and {L'^*) 
link configurations, the two-point function is perfectly 
real and the probability density is proportional to the 
square of the wave function. In this analysis, we choose 
to look at the zero-momentum probability distributions 
3 time slices after the source. 

We use the 2 + 1 flavour 32^^ x 64 PACS-CS config- 
urations 1 22 1, constructed with the Iwasaki gauge ac- 
tion [26 J with p - 1.90, giving a lattice spacing of 
0.0907(13) fm, and the (9(a)-improved Wilson action 
||27 l. We use 198 gauge field configurations, and em- 
ploy multiple sources per configuration, separated by at 
least one quarter of the temporal lattice extent. The hop- 
ping parameter for the light quarks is K^d = 0.13781, 
giving a pion mass of 156 MeV. 

To cleanly access the first three states, a 4 x 4 varia- 
tional basis is constructed using the;^'i operator with 16, 
35, 100 and 200 sweeps of Gaussian smearing [281, cor- 
responding to RMS smearing radii of 2.37, 3.50, 5.92 
and 8.55 lattice units respectively. We fix to Landau 
gauge by maximizing the 0{a^) improved fixing func- 
tional 1^ 



Tnnp = y Retr(-(f/^W- -— (f/^(x)t/(x + A) + h.c.) 

^—' i iZUn 

(10) 
using a Fourier transform accelerated algorithm fi3U\ . 

The wave functions observed for all our states show 
an approximate symmetry over the eight octants sur- 
rounding the origin. To improve our statistics we av- 
erage over these eight octants before presenting the re- 
sults. 

Our point of comparison with previous models of 
quark probability distributions comes from a non- 
relativistic constituent quark model with a one-gluon- 
exchange motivated Coulomb 4- ramp potential. The 
spin dependence of the model is given in Ref. |3 1 and 
the radial Schrodinger equation is solved with boundary 
conditions relevant to the lattice data; i.e. the derivative 
of the wave function is set to vanish at a distance LJ2. 

The ground state probability density for the d quark 
about the two u quarks at the origin obtained in our lat- 
tice calculations is illustrated in Fig. [2] We see that 




Figure 2: The probability distribution of the ground state d quark 
about the two u quarks fixed at the origin. 



the well-known sharp-peaked shape associated with the 
Coulomb potential is reproduced. 

The lattice data are compared with the constituent 
quark model in Fig. [3] Here, both the quark model 
probability distributions and the lattice results have been 
scaled such that the peak value is 1. The two quark 
model parameters adjusted in the fit are the string ten- 
sion, Vo" = 440 ± 40 MeV, and the constituent quark 
mass, niq ~ 370 MeV (accommodating the fact our 
quark mass is above the physical value). Using a least- 
squares fit varying the parameters nig and y/cr, we find 
the ground state lattice results are described well with 
s/cr - 400 MeV and niq - 360 MeV, which gives a 
ground state mass of 940 MeV and a first excited state 
mass of 1573 MeV. These parameters are held fixed in 
examinations of the excited states. 

Lattice results for the d-quaik probability distribu- 
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Figure 3: Comparison of the ground state lattice probability distribu- 
tion with the quark model. The probability distiibutions are qualita- 
tively similar, as would be expected for the ground state. 




Figure 4: The probability distribution of the d quark about the two u 
quarks at the origin in the first excited state. The darkened ring around 
the peak indicates a node in the probability distribution, consistent 
with a 25 state. 




Figure 5: The isovolume of the probability distribution of the d quark 
in the first excited state (colour map as in Fig.|4l. The outer edge can 
be seen to be affected by the boundary, indicating a necessary finite- 
volume effect associated with multi-particle components of the state. 



tion and the finite volume of the lattice. At this very 
light quark mass the distortion of the probability den- 
sity is significant and will correspondingly influence the 
lattice hadron mass. This interplay is key to extracting 
resonance parameters from lattice simulation results. 

Comparing the lattice probability distribution for the 
d quark in the first excited state to that predicted by 
the constituent quark model in Fig. l6j we see a quali- 
tative similarity but with important differences. While 
the node position is similar, the shape of the wave func- 
tion tail is different. 

The probability distribution of the second excited 
state of the nucleon in Fig. l7] reveals two nodes con- 
sistent with a 3S radial excitation of the d quark. Finite 
volume effects become even more apparent as shown in 
Fig. [8] distorting the outermost shell of the wave func- 
tion into an almost square shape. Comparing this state 
to the quark model prediction in Fig.l9] we observe qual- 
itative agreement. 

In this world-first study of the quark probability dis- 
tribution within excited states of the nucleon, we have 
shown that both the Roper and the second excited state 
display the node structure associated with radial ex- 
citations of the quarks. On comparing these proba- 
bility distributions to those predicted by a constituent 
quark model, we find good qualitative similarity with 
interesting differences. The discovery of a node struc- 
ture provides a deep understanding of the success of 
the smeared-source/sink correlation matrix methods of 
Ref. 113. 

Finite volume effects were shown to be particularly 
significant for the excited states explored herein at very 
light quark mass. As these excited states have a multi- 



tion in the first excited state of the proton are presented 
in Fig. p] The distribution exhibits a hydrogenic node 
structure consistent with a 2S state, indicating that the 
state includes a radial excitation of the d quark. This 
structure also indicates that the ideal combination of op- 
erators to access this state on the lattice would be super- 
posed Gaussians of different widths and opposite signs. 
This observation validates the approach of combining 
multiple smearing levels to construct the variational ba- 
sis and indeed the alternating signs of superposed Gaus- 
sians are observed in Refs. ifTSllSTI . 

The isovolume of this probability distribution illus- 
trated in Fig.lSlclearly shows the nodal structure, with an 
inner sphere surrounded by a near-spherical shell. The 
deviation from spherical symmetry in the outer shell di- 
rectly displays the important interplay between the en- 
ergy of the excited state observed in the lattice simula- 
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Figure 6: Comparison of the first excited state rf-quark probability 
distribution from our lattice QCD calculation (crosses) with the quark 
model (solid curve). The quark model predicts the node in approxi- 
mately the correct location, but deviates at the boundary. 




particle component, the interplay between the lattice 
volume, the wave function and the associated energy are 
key to extracting the resonance parameters of the Roper. 
Future calculations will explore the structure of the 
Roper in more detail, examining the mass dependence 
of the wave functions, more general spatial configu- 
rations of the quark positions, and the introduction of 
isospin-1/2 spin- 3/2 interpolating fields to reveal the 
role of D-wave contributions to the Roper. 
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Figure 7: The probability distribution of the d quark in the second 
excited state of the nucleon. Two nodes are visible, consistent with a 
35 state. 




Figure 8: The isovolume of the probability distribution of the d quark 
in the second excited state. The outermost node is compressed by the 
boundary into an almost square shape, indicating strong finite-volume 
eifects. 
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Figure 9: Comparison of the second excited state d-quark probabil- 
ity distribution from the lattice (crosses) with the quark model (solid 
curve). The nodes in the lattice data fall in between those predicted 
by the quark model. 
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